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Semi-analytic treatments of the evolution of orbits of weakly interacting massive particles 
(WIMPs) in the solar system suggest that the WIMPs bound to the solar system may enhance 
the direct detection rate relative to that of the unbound population by up to a factor of order 
unity, and boost the flux of neutrinos from WIMP annihilation in the Earth by up to two orders 
of magnitude. To test these important but uncertain results, we perform a suite of numerical orbit 
integrations to explore the properties of the bound WIMP population as a function of the WIMP 
mass and the scattering cross section with baryonic matter. For regions of WIMP parameter space 
presently allowed by experiments, we find that (i) the bound WIMP population enhances the direct 
detection rate by at most ~ 1% relative to the rate from unbound halo WIMPs; (ii) it is unlikely 
that planned km 3 -scale neutrino telescopes will detect neutrinos from WIMP annihilation in the 
Earth; (iii) the event rate from neutrinos produced by WIMP annihilation in the Sun may be much 
smaller than implied by the usual calculations, which assume that WIMPs scattered onto bound 
orbits are rapidly thermalized in the Sun. 
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1. Introduction 

WIMPs in the solar system may be detected in (i) "direct detection" experiments, which mea- 
sure the recoil of nuclei during interactions with astrophysical WIMPs [jl|] and have placed inter- 
esting constraints on the WIMP mass m x and spin-independent (or -dependent) elastic scattering 
cross section with protons Op 1 and neutrons o^ 1 (<Jp D and C„ D ); and (ii) "indirect detection" exper- 
iments, including neutrino telescopes. The next generation of neutrino telescopes may detect the 
high-energy neutrinos from annihilation of WIMPs captured in the Sun and Earth [Q]. 

For a given WIMP model, event rates in direct and indirect detection experiments are deter- 
mined by the phase-space distribution function (DF) of WIMPs. The fiducial assumption is that 
these event rates are dominated by WIMPs from the Galactic halo, passing through the solar sys- 
tem on unbound orbits [ph. However, the following two processes may create additional long-lived 
populations of WIMPs bound to the solar system, (i) Damour & Krauss [Q] argued that secular 
gravitational interactions with the planets could increase the perihelia of Galactic WIMPs scattered 
in the Sun to larger radii, quenching further scattering by solar nuclei and dramatically increasing 
the WIMP lifetimes. Their semi-analytic estimates indicated that the bound population could in- 
crease direct detection rates by up to a factor of ^(1), and enhance the neutrino flux from the Earth 
by two orders of magnitude [Q]. (ii) Gould and Lundberg & Edsjo ^ estimated that the density of 
bound WIMPs at Earth due to gravitational capture and scattering by the planets was ~ 1 % of the 
density of unbound WIMPs if additional elastic scattering in the Sun were unimportant. Because 
of their low geocentric speeds the bound WIMPs could enhance the capture rate in the Earth by up 
to an order of magnitude, hence boosting the neutrino flux. 

Semi-analytic techniques cannot, however, describe the full spectrum of behavior in few-body 
gravitational systems like the solar system. Numerical orbit integrations are the only reliable tool 
for this task. In this paper, we summarize the results of a set of such integrations, which are 
described further in 

2. Simulations 

We assume that the DF of the halo dark matter is Maxwellian, with a local halo WIMP number 
density n x = p x /m x . m x is the WIMP mass, and we set the local WIMP mass density p x = 
0.3 GeV cm~ 3 . We set the speed of the Sun relative to the halo to v = 220 km s , and the one- 
dimensional WIMP velocity dispersion a = v Q / yjl. 

The initial conditions for the orbit integration are of two kinds. Type I: We employ Monte 
Carlo methods to determine the positions and velocities of bound WIMPs after they have scattered 
elastically off nuclei in the Sun (using the solar model in 0). We run four sets of simulations, each 
corresponding to a particular choice of WIMP spin-independent elastic scattering cross section and 
mass {qf = 10~ 41 cm 2 : m x = 60 AMU ; af= 10~ 43 cm 2 : m x = 60, 150 and 500 AMU). The 
spin-dependent cross section is set to zero for simplicity. Type II: To estimate the DF of WIMPs 
captured gravitationally, WIMP orbits are followed numerically after they approach within 1000 
AU of the Sun. We carry out the simulations with only a single mass and cross section (m x = 500 
AMU, Op I= 10~ 43 cm 2 ), but record the integrated optical depth from passages through the Sun as 
a function of time, which allows us to extend the results to other masses and cross sections. 



2 



WIMPs in the solar system 



Annika H. G. Peter 



Geocentric: = 500 AMU, of = 1 43 cm 2 



Unbound 

Spin-Independent Capture ■ — ©■ 




v [km/s] 

Figure 1: DFs as a function of geocentric speed v, normalized so that / f(v)v 2 dv is the number density; n x 
is the halo number density in the solar neighborhood. The dashed line is the unbound halo DF. The circles 
represent Type I initial conditions with <Jp D = 0, the squares show the bound DF from gravitational capture 
(Type II), and the solid line without points is the estimated DF for Type I capture with <Jp D = 10~ 36 cm 2 . 

Bound WIMPs are followed numerically, admitting the possibility of additional elastic scat- 
tering in the Sun. We use a simplified solar system containing only one planet, Jupiter, revolving 
on a circular orbit about the Sun. The integration algorithm is symplectic with an adaptive timestep 
determined by the heliocentric distance We use other methods to integrate orbits when WIMPs 
go through the Sun or have close encounters with planets. Integration is terminated before the age 
of the solar system if the WIMP escapes from the solar system or is scattered in the Sun onto an 
orbit with semi-major axis a < 0.5 AU. Type I simulations followed 10 s — 10 6 particles. Type II 
simulations followed over 10 10 particles, of which ~ 3 x 10 5 were at least temporarily bound. 

3. Results 

We present a sample of our results in Fig. [I], which shows DFs f(v) at the Earth as a function 
of geocentric speed v, in the case of m x = 500 AMU and <Jp I= 10~ 43 cm 2 . The dashed line denotes 
WIMPs from the halo, unbound to the solar system. Blue circles show the DF of WIMPs with Type 
I initial conditions with Gp D = 0. Red squares represent the WIMPs with Type II initial conditions. 
Errors are estimated using bootstrap resampling. We extrapolate our results to obtain a maximum 
DF due to spin-dependent scattering (the solid magenta line). The main conclusion in Fig. [j] is that 
the total bound number density is small, no more than ~ 0. 1% of the unbound density. 

We estimate the direct detection rates of bound WIMPs, assuming that spin-independent inter- 
actions dominate in the detector volume. The analysis thresholds of most experiments lie above the 
maximum possible recoil energy for bound WIMPs (< 10 keV). An exception is the XENON10 
experiment [jl|], for which the maximum enhancement to the differential event rate is ~ 1% of the 
halo event rate. The enhancement to the total event rate in the analysis window is ~ 0. 1%. 
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We determine event rates in neutrino telescopes from WIMP annihilation in the Earth using 
the DarkSUSY software package, assuming the WIMP is a supersymmetric We find that the 
event rate in an IceCube-sized neutrino telescope [§] is below the detection threshold. 

The Sun is a target for WIMP annihilation searches. Initial results from the Super- Kamiokande 
experiment suggest that neutrino telescopes will be competitive with direct detection experiments 
in their sensitivity to Op D , with the standard assumption of rapid thermalization of WIMPs in the 
Sun JTo|]. The neutrino flux from the Sun is usually estimated assuming that all particles scattered in 
the Sun onto bound orbits rapidly reach thermal equilibrium with solar nuclei, settling into a dense 



core at the center of the Sun if m x > 4 Gev [11]. Thermalization is crucial since the annihilation 
rate goes as the square of the density. This picture may be realistic if m x < 1 TeV. However, we 
find that most WIMPs on Jupiter-crossing orbits are ejected from the solar system before they can 
rescatter in the Sun (so long as Op D < 10~ 38 cm 2 ). There are also long-lived WIMP populations 
that may not thermalize. We find that if o s p D < 1(T 43 cm 2 (> 1(T 43 cm 2 ) and m x = 1 TeV, 20% 
(12%) of WIMPs will never thermalize in the Sun, and rises to near 100% (85%) for m % = 10 TeV. 
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